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INTRODUCTION 


This report is concerned with the application of the elements of 
quantum (wave) mechanics to some special problems in the field of macro- 
scopic fluid dynamics (mechanics), often referred to as classical fluid 
dynamics. In particular, considerations will be on flow of a viscous 
imcompr es sible fluid around a circular cylinder. The presentation is 
divided into three sections: Section 1 constitutes a brief presentation 
of the flow of a nonviscous fluid around a circular cylinder. Attention 
is called to the fundamental concepts in any frictionless fluid flow such 
as velocity potential, stream function and so on. Section Z presents a 
brief discussion of the restrictions imposed upon the stream function by 
the number of dimensions of space in which one usually operates (two in 
contrast to three), by the differences between stream function in two- 
and three-dimensional flows, and by the selection, as in the present 
case, of the two-dimensional space representation. The third section 
presents in detail the flow past three-dimensional bodies in a viscous 
fluid, particularly past a circular cylinder in the symmetrical case. 
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1. FLOW AROUND THREE-DIMENSIONAL BODIES--FRICTIONLESS 
FLUID FLOW 

1.1. Characteristic Functions of the Flow 

The part of fluid flow theory dealing with the characteristic functions 
of the non-viscid frictionless fluid flow with no heat conductivity is 
usually well understood by most scientists, physicists, engineers, and 
chemists and very little remains to be added. It appears desirable to 
call attention to a few of the most characteristic laws and functions which 
are valid and widely used in the field of frictionless fluid flows. Subse- 
quently, these laws and functions are extended to and used in the field 
of viscous and heat-conductive fluid flows with possibly small modifications 
and adjustments, wherever they may become necessary. The most 
significant laws, rules and functions in the frictionless fluid flows are: 

a. conservation of mass (or the continuity equation); 

b. conservation of momentum; 

c. irrotationality; 

d. velocity potential; 

e. relation between the irrotationality and the existence 
of the potential function; 

f. conservation of energy. 

The stream function requires a different definition for two-dimensional 
flows and a different one for three-dimensional flows. A streamline, 
however, is defined in the same manner for either two- or three-dimensional 
flow; viz. , a continuous line through the fluid such that it has the direction 
of the velocity at every point throughout its length. The differential 
equation for a streamline in three-dimensional flow is: 

dx/u = dy/v = dz/w; (1. 1. 1) 

whereas in two-dimensional flow the form is: 

dx/u = dy/v, (1. 1.2) 

In two-dimensional flow the concept of the stream function, 4^ =4^ (x,y), gives 

u = -^/8y; V = 844/9 x; 9u/9y = 9v/9x; (1. 1.3) 
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(1.1.4) 


and finally for the irrotational flow cases 

+ 8^1^ /8y^ = V^\|^ = 0, 

2 

where the symbol v denotes the Laplace operator (Laplacian). 

2. SOME CHARACTERISTIC FUNCTIONS FOR TWO- AND THREE- 
DIMENSIONAL FLOWS 

2. 1. Stream Function in Two-Dimensional Flow 

The stream function requires a different definition for two- 
dimensional flow from that for three-dimensional flow. One of the most 
useful definitions of the two-dimensional stream functions is that the 
partial derivative of the stream function with respect to any direction is 
the velocity component plus 90 degrees (counterclockwise) to the direction 
of flow (Streeter, p, 39). In addition it must satisfy the continuity equation. 
The first part of Equation (1.1,3), i, e. , u = . . . , v = . . . , is true whether 
the flow has rotation or not. For irrotational flow , however, the second 
part of Equation (1. 1.3) and Equation (1,1.4) are true. In particular, 
Equation (1. 1.4) shows that the stream function may be constructed as 
the velocity potential for some other flow (Streeter, p. 39). Relations 
between stream function, ijj, and velocity potential, 4>, are found by 
equating the expressions for the velocity components: 

8<|)/8x = 8i|i/8y; 8c|>/8y = -d^j/Ox. (2. 1, 1) 

2. 2. Stream Function in Three-Dimensional Flow 

Again, all of the details of three-dimensional flow will not be 
presented, only those which may be pertinent to the specific problem of 
the present report. The calculation of the stream function in three- 
dimensional flow requires much more effort and energy than that for the 
two-dimensional case. Of special interest and value is the Stokes* 
stream function defined only for those three-dimensional flow cases which 
have axial symmetry. 

The above discussion demonstrates clearly the advantages of using 
the two-dimensional flow representation over three-dimensional flow 
representation, even in the cases where the actual flow in the physical 
space is geometrically expressed as a space phenomenon. The following 
representation of the flow phenomena in the viscous, incompressible 
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fluid in a three-dimensional space configuration will only use those 
analytical methods which give the mathematical description in terms of 
two independent variables. 

3. FLOW PAST THREE-DIMENSIONAL BODIES IN VISCOUS 
INCOMPRESSIBLE FLUID 

3. 1. Flow Past a Circular Cylinder --Symmetrical Case 

Flows past three-dimensional symmetrical bodies, in particular 
flow past a circular cylinder in a symmetrical configuration, will be the 
first case to be considered. The flow in the boundary layer past a 
circular cylinder in the symmetrical case is known from the literature, 
and the solution functions used here are taken from Schlichting, pp. 150-152 
The free stream velocity is parallel to x-axis , and one begins with the 
ideal velocity distribution in "potential" irrotational flow past a circular 
cylinder of radius R and free stream velocity U^. In general, when one 
considers the flow past a cylinder in the symmetrical case, one may use 
the so-called Blasius series. According to Blasius. the potentiaT flow 
is given by the series: 

3 5 7 9 ll 

U(x) = u^x + u^x + UgX + u^x + u^x + Uj^x + . . . ; (3. 1. 1) 

where the coefficients u^, depend only upon the shape of the body 

and are to be considered known. The pressure term in the boundary layer 
can be easily calculated for the stationary condition and becomes: 

- dp/dx = U dU/dx ; p = constant; (3. 1.2) 

2 3 5 2 

U dU/dx = UjX + 4u^u^x + x (6 u^^u^ + 3 u^) 

+ x^(8 Uj^u.^ + 8 u^u^) + x^ (10 + 10 ^3^.^ ^ ^^5^ 

+ (12 +12 ^3^^^ +13 + . • • ; (3. 1. 3) 

where the continuity equation is integrated with the use of the stream 
function ijj = i|i(x,y). A suitable assumption for the coefficients of the 
stream function and hence for the velocity components remains to be 
made. For the general case, Blasius proposed the form for k|;(x,y), u 
and V with the standard boundary condition (y=0, u = v = 0, y=“. 
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The case treated in this work refers to the boundary layer on a 
cylindrical body placed in a stream which is perpendicular to its axis 
(this axis will be denoted as the z=axis, perpendicular to the plane (x,y) 
in which the flow actually takes place). Thus the flow is considered to 
be two-dimensional in the (x,y) plane. In either case the velocity of the 
potential flow is assumed to have the form of a power series in x, where 
'‘x” denotes the distance from the stagnation point measured along the 
contour of the cylinder. In the Blasius series, the velocity profile in 
the boundary layer is also represented as a power-series in x, with the 
coefficients assumed to be functions of the coordinate y, measured at 
right angles to the wall; i. e. , to the surface of the cylinder. 

The form is ; 

U(x) = 2 sin 4> = 2 sin (x R ^), (3. 1.4) 

and by expanding in a power series it becomes: 

U(x) = 2 [ (x R"^) » i (X R“V + i (x R' V 


- i (x R ^)^ + . . . ] ; 

the stream function with r| = y R ^ (2 U^R ^)^^ is: 

4i = (i^u“^) {u^x f^(ri ) + 4 u^x^ f3(ri ) + 6 u^x^ f^ (ti ) 

+ 8 u^x^ f.^(Ti ) + 10 u^x^ fj^(ri ) + 12 f^^(-n ) + . . . } 


(3. 1.5) 


the functions u are: 
n 

1 ^323 

^1 ~ ^ * ^3 “ “ Tl ~ *51 ; . . . ; 

and the velocity u is given by: 

u = 2 {(xR‘^ f'j - (i ) (xR'V f'3 + (|, ) (xR‘S^ £'5 


{j, ) (xR'^”^ + . . . } . 


where the functions f are called the functional coefficients, 

n 


6 



In order to render the functional coefficients independent of the 
particular properties of the profile, i. e. , of Uj^, u^, . . . # it is necessary 
to split them as follows ( for i « 5): 

f'j, = g'5 + (10/3) h'^; (3.1.9) 

= g'^ + 7 + (70/3) k-^; (3.1.10) 

f'g = g'g + 12h'g + (126/5) k'g + 84 j'g + 280 q'g: (3.1.11) 

f'll = g'jj + (155/3) + 66 k'jj + 220 j'jj + 462 q’jj . . . ; (3. 1. 12) 

where all the functions shown are tabulated in Schlichting's work. The 
usefulness of the Blasius' method is restricted by the fact that it should 
not be used for slender bodies. Other attempts on the problem have been 
made by Goldstein, Hiemanz, Goertler and others. For a circular 

cylinder of diameter d = 9. 74 cm in a stream of velocity u = 19 cm sec , 

-14 ^ 

Re = u d 1 / = 1. 85 X 10 , Hiemanz found that his experiment gave values 

o 

for u which could be represented sufficiently accurately by three terms 
of the series for the potential flow. Due to the fact that the function f^(r| ) 
are tabulated in Schlichting's work, this research will use the Schlichting 
data throughout the report (see Schlichting, p. 152 to p. 154 after remodeling) 

3. 2. Forms Used in the Analysis 
Stream function 

(2 V {(R'^x)fj(ri) - |j (R'^x)^ f 3 (r,) 

+ (R- 1 x)^ in) - fy (R‘ ^ + -^ (R‘ ' X)'’ fg(n) 

- (R“ (3. 2. 1) 

Velocity Components (Schlichting Data) 

The velocity components are derived from the stream function: 

u = (94j/8y) = (84>/9n) • <2^^^ v 
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= 2 {(xR'^f'j(T7) - ^ +|l (xR-^^f^tr,) 

- It (xr' + If (xR‘ i^(T)) 

- (xR" + « o , } c (3o 2. 2) 

Similarly, the velocity component v is: 

V = - 84j/8x ; 


V = - (2 V U^R)*/^ {R'^fj(r)) - |t R‘^x^ £3(11) 

+ -Ij R"®x^ £5(11) - fi R'^x^ {y(n) + If 

- ■— R‘ £jj(ti) +...}; 71 = (2 R" ^ 1 /" y, 

(3. 2.4) 

In Schlichting's work one can find another form for the function v, (p. 149, 
Equation (9 o 20a)) which, after inserting the corresponding values for the 
coefficients u^, is in perfect agreement with Equation (3.2,4). 

Final Forms of the Vorticitv Component 


3v/3x = - (2 V U„ R)^/^ ^ If 

+ (.^)(xR-l)^ f (n) - (^)(xR'h^ fj(n) 

+ (■^)(xR‘^)’^£9(n) - (-^f^)(xR'‘)'^fii(n) + ••• } 

(3. 2. 5) 


Ou/&y) = 2 U^(2 R'^ {(xR'^> £'j'(77) 

- (f-, ) (xR‘ fy(T)) + (-|r) (X R' E^'(r) 

- e^) (xR‘^)^fy(ri) + (|f) (xR'*)^ fgMti) 

- (-^)(xR‘^)“ £'j’j(ri) + ... } • 


(3. 2. 6) 
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The function o) in the form 
z 


w =i [{0v/3x) - 8u/9y)] , » sec , 

Zi u 


(3. 2.7) 


-1 


will be superimposed upon each curve u to illustrate the perturbed 

velocity distribution in the laminar boundary layer on a circular cylinder 
at the following set of angles: 

<!> = 0°; = 20°; = 40°; = 60°; = 80°; = 90°; = 100°; = 110°. (3.2.8) 

where = (x R ^). 


4. PERTURBATIONS 
4. 1, Stream Function 

Thus, there have been derived the equations for the stream function, 
= k|;(x,Ti ) = 4^(x,q{y)), and the equations for the velocity components, 

(u,v), in the specific flow in question, i.e. , the symmetrical case of 
the flow past a circular cylinder (Blasius series). From the definition 
of the stream function one can derive the equation of the streamlines. 
Since, the stream function satisfies the equation of continuity, the curves, 
= constant, represent streamlines (Owczarek, p. 63). The research 
begins with the stream function (see Equations (3.2. 1) and (3,2.4) above): 

i), = (2 1/ {(R'^x)fj(T, ) - i (R'^x)^ f3(ri ) 

+ (R'^x)^ - j, (R‘V f7(9) +-^ (R‘ x)^ 

_ (R”^x)^ ^ f^ (r| ) + .... } = constant. (4.1.1) 


Put in this form, the function i|i is a function of two variables, (x,r|), or 

(x, T]),x = x r”^ (dimensionless), or (x,y), since r| = constant y, where 

the constant is a dimensional number. Thus, both the coordinates 

(x, q) appearing in the representation of the function i|j, Equation (4. 1. 1) 

are dimensionless , and only the coefficient (2 i/ U standing in 

2 -I . 

front of { . . . } has the dimensions of cm sec as it should be. The 
part of the stream function in the braces, { . . . } , must be kept constant 
in order that the stream function be constant (condition for the streamlines). 
In other words, the condition for streamlines appears in the form: 


9 



(4, 1,2) 


(R"^x) (n) - 7 ; (R'M^ fjC’l) + I; ^<’1* 

,, i (R^'^x)^ f^(r| ) + . , , , = constant; 
and one may group the above terms in the following manner: 

[ (R°^x) (ti )] ■=■ [ j, (R“^x)^ £3(11 )] +[ji fgh)] 

. [ (r‘V f7(>i)] + (R‘M^ ^ 9(^)3 * [tt: 

= constant = C , 

Obviously, there are infinitely many possibilities for dividing the constant 
appearing on the right hand side of Equation (4. 1. 3) among the infinitely 
many terms of the series appearing on the left hand side of the same 
equation. For illustrative purposes, the following scheme is proposed: 
the constant C is divided into "m" equal parts, where "m" is equal to 
the number of terms on the left hand side of (4.1 . 3); each term is assumed 
to be equal to C m and (R ^ : 

5fj(ti) = C m'^ 

- y, (x)^ f j(tl ) = C m S 

(x)5 £ 5 ( 9 ) = C m‘S (4.1.4) 

-i (J)''f7(ri) =Cm-S 

) = C m j 

and so on 

although there may be other alternatives, one of the conclusions of the 
system of equations, (4. 1, 4), is that, in the first form of (4. 1, 4), both 
the variables, (x,r|), or any logical combination of these variables such 
as the product or the ratio are assumed to be constant. Again from the 
many possible combinations, which may and probably do occur in the 
actual, physical conditions, the following combination are selected for 
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purely illustrative purposes: 


X Ti = constant = C^, 1.5) 

[ ©".1 ] = constant; t; = y R"^ (2 U„Ri-'V^^ (4- 1- &) 

The first possibility. Equation (4. 1. 5), will be discussed in more detail. 

As can be seen from Equation (4. 1. 5) this is an equation of a hyperbola. 
Consequently, in the case under consideration, the simplest possible 
geometrical form of a streamline in the (x,y) space is a hyperbola. 

Equation (4. 1, 5). 

4.2. Elementary Geometrical Characteristics of the Hyperbola 

The geometrical characteristic of hyperbola are collected and presented 
in oraer to use tnem in Uie present work, thus eliminating the need for 
outside references. The standard form of a hyperbola is: 



where the axis, x, intersects the hyperbola at the vertices and 
the segment A A^ is the transverse axis; both branches of hyperbola 
are symmetrical with respect to (A^ A^), one part being below and another 
above the axis; the y axis does not intersect the hyperbola since the center 
of the coordinates (x,y) is located in the center between A^ and A^. 

From Equation (4, 2. 1), if a = b, the hyperbola is equilateral, and if the 
asymptotes are perpendicular to each other, the hyperbola is "rectangular. " 
The equation of a rectangular hyperbola in which its asymptotes are 
referred to as coordinate axes is: 

x y = a^ = A , (4. 2. 2) 

= 2A . 


and is located in the first and third quandrants. For such a hyperbola 

12 1 2 
one has x y = a (first and third quadrants), (second and 

fourth quandrants); 


semi-axes : 


a = b = (2A)^/^ ; 

c^ = + b^ = 4A; c = 2A^/^; 


(4. 2. 3) 
(4. 2.4) 
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J/2 


. . -1 , 1/2 

eccentricity I e = c a - i 


1/2 


1/2 


0) : 

(4. 2. 5) 

: [(2A)^/^. 0] : 

(4. 2. 6) 


(4. 2.7) 


(4.2.8) 

X = .a1/=^ ; 

(4. 2. 9) 

X + y = 0; or X = y, X = 

-y (4.2.10) 


Consequently, equation (4. 1.5) is the equation of a rectangular hyperbola 
in which its asymptotes are referred to as dimensionless coordinate 
axes (R"^x,r|) or (x, q ). This coordinate system will be used as the 
basic coordinate system in further considerations, discussions and 
plotting of the diagrams. 


4. 3. Computer Plots 

The plots, included in the present report refer to a solution of the 

steady-state boundary layer equation in the flow past a circular cylinder 

in the symmetrical case (Blasius series). In particular the plots refer 

-1 

to the horizontal velocity component in the case given by u as quoted 
previously? 


u = 2 { X f'^Cq ) “ ~ (x)^ f'^ (q ) + . . . } (4.3.1) 

-1 

The plots of the function of u as a function of two variables are given 

in the Schlichting book (p. 153) where the function u is a function of 

y (U^ Ri/”^)^^^ = (q 2“^^^) and ofcj) = xR"*^, with ^ being treated as a 

parameter in the range of = 0°, 20°, 40°, 60°, 80°, 90 , 100 , 108. 8 , 

Similar plots have been accomplished from this research which are done 

automatically in their entirety from the output plotter as part of the 

CDC 650 Computer System. These plots are functions of the independent 

variable q varying from q = 0 up to q = 4, 0, and for the values of the 

-1 

parameter c|), (j> =xR : 


(}) = 10°. 20°. 30°. 40 , 50 . 60 , 70°. 80°. 90 . 100 


no 


Since this effort deals with the problem of the flow in a viscous 
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fluid (viscosity is defined as a transverse transfer of momentum) any 

possibility of obtaining a purely theoretical laminar flow around a 

cylinder is automatically excluded. One can only talk about a quasi- 

linear flow or a flow with small disturbances due to small viscosity 

pheomena. In order to make the disturbances due to viscosity visibly 

more observable, the concept of disturbance in the form of the vorticity 

function, w , containing (8v/8x), and (8u/8y), has been introduced which 
z -1 . . 

is then superimposed upon the velocity function, u . This is done 

geometrically, by means of plotting the diagrams of (u ) as the 

function of q in the two-dimensional Cartesian coordinate system, 

(u and q = y (2 R with (f) = x = (x R~^ being a 

parameter, and R being a fixed constant. Moreover, the function 

is also plotted in the same diagram and in the same coordinate system 

as a function of r| . Under the assumption that one is dealing with the 

linear Schroedinger equation as the basic equation, the principle of 

superposition can be applied and the values of the two functions u U^^^ 

and (j , and be added geometrically for each value of the independent 

variable, q . The final operation can then be performed for obtaining 

the sum of functions. This technique has been applied and the resulting 

plots are included as Appendix A of the report. It is not essential which 

coordinates in all of the systems discussed be chosen: (x,y) or (x,y), 

(x, T|), (x,r|), since obviously the coordinates x and x, y and q differ 

only by constants, one being dimensional and the other dimensionless , 

respectively. In conclusion, one may summarize the results obtained 

as follows ; 

(a) the stream function of the flow around a circular cylinder, 
located symmetrically, is given by Equation (4, 1.1); 

(b) the velocity component u is given by Equation (3. 2. 2); 

(c) the velocity component v is given by Equation (3. 2, 4); 

(d) the disturbance function (in form of the curl of velocity) is 
given by Equations (3, 2, 5), (3. 2. 6), and (3. 2, 7); this is 
(geometrically) superimposed upon the velocity component 
function u U^^; 

(e) the condition for the streamlines is given by Equation (4. 1. 3); 

(f) the simplest possible geometrical form of a streamline in 
the (x,y) or (x, q ) space is the hyperbola (xq ) = C^; 
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(g) consequently, the result of the case of flow around a circular 
cylinder is exactly the same as in the previous case of the 
laminar flow along an infinitely long plate. A finite but very 
small amount of vorticity introduced into the flow system 
causes the appearance of small disturbances. These disturbances 
originate at each finite element and spread in the fluid for short 
distances. In the case of the flat plate the streamlines were 
parabolas whereas in the present case they are hyperbolas 
as proven in Equation (4» 1. 5). 

As in the case of a flat plate, one assumes a moving coordinate 
system (x,y) or rather [ (xR"^),q] along each of the curves given as the 
sum of [ (u U”^) + W ] , from Equation (3. 2. 2), plus Equation (3. 2. 7) with 
Equations (3. 2. 5) and (3. 2. 6), supplying a series of values of the para- 
meter =0®, 20^, . . . , 110®. Each chosen point for a certain 

value of r\ is assumed to be the vertex of the rectangular hyperbola. 

The moving coordinate system moves parallel to itself and to the fixed 
coordinate system, [ {xR”^), vertical axis] . The branches of the hyperbola 
are located in the first and third guadrants. From both coordinate axes, 
the horizontal axis (R x) is always horizontal and the other axis is 
always vertical. By means of the proper transformation (of coordinates 
if necessary) the branch of hyperbola located in the third quadrant is 
being shifted along the transverse axis (Aj^ - A^) so that both vertexes 

A and A coincide at the vertex of the branch of the hyperbola located 
X 2 

in the first quadrant. Consequently, the first quadrant branch of hyperbola 
is oriented from up-down and to the right; the third quadrant branch of 
hyperbola is oriented from down— up and to the left; both branches meet 
at their vertexes (A^ at the same point as A^) and both vertexes coincide 
with the point chosen on the curve [ (u + u)^]. The two branched hyperbolas 

in question are traced at a number of points along the curves [ (u 
and their points of intersections are located. The sections of streamlines 
(hyperbolas) between two found points of intersection furnish the zig-zag 
pattern of the path of a particle. This produces the resulting pattern of 
disturbances in the flow around a circular cylinder caused by the injection 
of the vorticity geometry into the regular pattern. These zig-zag patterns 
are traced in a few cases. The data used in the computer plots are the same 
referred to previously; 
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U = 200 km hour ^ = 5555. 55 cm sec S 

00 

Z -1 
1 / = 0. 149 cm sec . 

4. 4. The Plots 

The first set of plots refers to the flow around a circular cylinder 
in the symmetrical case. The velocity distribution, (u ), is taken 
as a function of (x for various values of the coordinate r\ , considered 

to be a parameter. The value of t] for 20 values is taken to be (t] - 

Ti =-0. 2, 0.4, 0. 6, 3. 8 = 4. 0; (4.4. 1) 

In this approach no disturbances are superimposed upon the flow. The 
20 plots included at the end of the present report are based on the values 
of the function of (u U’^) as taken from Equation (3. 1.8) or (4. 3. 1). 

Por a certain number of points of these functions the streamlines are 
plotted as double-branched hyperbolas , and the points of intersections 
of these hyperbolas are found; the sections of streamlines between the 
points of intersection are interconnected thereby displaying the zig-zag 
pattern as the path of a particle. The results show that the zig-zag 
pattern, which indicates the existence of some disturbances in the boundary 
layer, appears even in the flow without the introduction of any outside 
disturbance such as vorticity. This can be expected in this kind of flow. 

It seems desirable to point out that one may question the accuracy and 
the precision of the graphical results and of the plots in the neighborhood 
of X R"^ equal to 2. 0. This may be partly justified since it must be kept 
in mind that: (a) the present calculations were done primarily for 

illustrative purposes, and for providing general conclusions giving 
insight into the particular points and pecularities of the problem needing 
more emphasis, and (b) the point in question (R ^ x = 2) lies outside the 
region of the separation of the flow and consequently the characteristic 
properties of the flow at this particular point cannot be measured with 
precision. 

The second set of plots is based on the curve u plus a disturbance 
with resultant curve; 

u + W = resultant curve; (4, 4, 2) 

CO Z 
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on each curve (u UjH taken from the first set of plots there is super- 
imposed the disturbance curve in the form of as the function of r\. 

Each graph therefore clearly shows two curves: one describing the value 
nf u as the function of and the second curve describing the vorticity 
function, u) , again as the function of t| . The vertical coordinates on 
these two curves are geometrically added giving the resultant curve, 

Equation (4. 4. 2). Again at a certain number of points on this curve, 
the two-branched hyperbolas are plotted and the points of intersections 
of the streamlines (i. e, , of hyperbolas) are located. The zig-zag 
patterns of the paths of particles are traced thus giving the final results. 

The third set of plots refers to the velocity v^ (dimensionless) = 

V Equation (3. 2. 4), which is subject to disturbances in the form of 

the vorticity, u , i. e. : 

V = V u"^; v from Equation (3. 2. 4); = resultant function 

^ (4.4.3) 

Again, at a certain number of points of the curve (v^ + u>^) two-branched 
hyperbolas are plotted, the points of intersections of the streamlines 
(i, e, , of hyperbolas) are located, and the zig-zag patterns of the paths 
of particles are traced, thus giving the final pattern. 

All the plots described above and included in this report demonstrate 
the existence of the zig-zag patterns in the flow in question, regardless 
of whether the flow is laminar, irrotational, rotational, or there is or 
is not a vorticity function geometrically superimposed upon the velocity 
functions of u and v components. The existence of the zig-zag pattern 
is proof that the flow is not laminar, the particle paths are not parallel 
(nor quasi-parallel) lines. This is an obvious indication that there are 
disturbances in the particular flow in question due to the shape (circular) 
of the body around which the flow takes place. The laminar flow along a 
flat plate demonstrates the zig-zag pattern. There does not exist a 
laminar flow in a domain of a viscous, heat-conducting fluid flow (above 
the \-transition point). This conclusion, which is in accordance with 
Heisenberg's statement, has been discussed in previous reports. The justi- 
fication for the use of summation law, association between the wave mechanics 
and deterministic macroscopic fluid dynamics, and all the other aspects are 
also explained in previous reports. 
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5. FINAL REMARKS 

A comparison of the geometry of the flow patterns of the laminar 
flow past a symmetrically located cylinder (Schlichting, pp. 146-155) and 
those in Technical Report No. 3 indicates a few important results: 

(1) The so-called "laminar*' flow around a cylinder {symmetri- 
cal case or not) is not a laminar but is a flow with disturbances; 

(2) Disturbances are always present in the so-called "laminar" 
flow due to the transverse transfer of momentum (viscosity); 

(3) The zig-zag paths appear always in the so-called "laminar" 
flow, however small they may be since streamlines are 
always there; 

(4) The "streamlines" in the "laminar" flow are actually the 
mean value paths of the real zig-zag paths of the particles; 
due to the fact that zig-zag paths are small (small amplitudes) 
and due to the physiological aspects they are seen by the 
naked human eye as continuous lines (the streamlines); 

(5) The above traced zig-zag paths seem to be very "regular" 
whereas the oscillograms obtained from the oscillographs 
placed in a turbulent jet (see Technical Reports No. 1 and 2) 
demonstrate that often the realistic zig-zag paths are 
irregular; a certain regularity in zig-zag paths may or may 
not appear in the "periodic" sense. This demonstrates that 
"local' irregularities (jumps, sharp steep "mountains", 

sharp, steep "valleys") have their origins and are due to 
other reasons and phenomena not yet discussed (interference, 
inter -correlation); 

(6) When the zig-zag paths are very "even", regular and of small 
amplitudes, a zig-zag path may be seen by the naked human 
eye as one thick "streamline", as one thick ray; a light, a 
ray of light, which is also a wave; 

(7) The problem of interference and inter -correlation is discussed 
in Technical Report No, 1. 
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The First Set of Plots 

The £\mction (u as the function of (x R for various 

values of the parameter r/ = 0, 2, 0. 4, . « . , 4. 0 . 

No disturbances. 
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The Second Set of Plots 


The function (u + u ) as the function of (x R“^) for various values 
' <D 2 

of the parameter r\ = 0. 2. 0.4 4, 0. The function u is 

dimensionless, the functions ^ sec Curl represents the dis- 

^ -*1 

turbances superimposed upon the function u . Since the dimensionless 
velocity and the disturbance have different dimensions, one has to assume 
that the plots are performed in time -frozen conditions. 
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The Third Set of Plots 

The function (v^ + disturbance), as the function of 

R for various values of the parameter t| = 0, 2, 0» 4, « . • j 
4.0, The function v , = v U is dimensionless, the disturbanc 
equal to u) i is superimposed upon the dimensionless velocity 
function, = v Time is frozen. 



ELOCITY DISTRIBUTION IN THE BOUNDARY LAYER ON A CIRCULAR CYLINDER 


KEY 

1 DISTUS^AHCE .bmOI 25IN 
I X/R a.RQI -- 3.5IN 
i UNIT HYPERBOLA - .LiCN . LblN 
ETA-^ 0.2 




DISTliRBANCF 



72 









74 



ELOCITY DISTRIBUTION IN THE BOUNDARY LAYER ON A CIRCULAR CYLINDER 

KEY 

i DISTURBANCE -- -- .25 IN 

i m B.ROi 

i UNIT HYPERBOLA -- .LICH ..IblN 
ETA-. 1.0 






CNi 







/ELGCITY DISTRIBUTION IN THE BOUNDARY LAYER ON A CIRCULAR CYLINDER 


KEY 

1 DISTURBANCE -- bCOT SIN 
i X/R - B.ROT -- 5.5IN 
i UNIT HfPEreaA MCW -- .TUN 
ETA-. 2.0 







DJSrUBANCt 



82 




83 









"VELOCITY DISTRIBUTION IN Ttf BOUNDARY LAYER ON A CIRCULAR CYLINDER 

KEY 


1 DISTURBANCE . .bUOl .251N 
I X/P, B.SCU 

i UNIT HYPERBOLA LiCM -- . IblN 
ETA. T.O 




MWA 




Q.u O.U Q.b 0.7 Q.B 0 

.0 I.O .i.L 1.2 I 

i.L i.O i. 

L. 1.7 l.B j.|q 

KANYA AND D. I*!IERENCA 


X/P. 

> 




DISTURBANCE 



86 













88 



1 


ELOCITY DISTRIBUTION IN THE BOUNDARY LAYER ON A CIRCULAR CYL 

KEY 

i DISTURBANCE -• .Id^aT -• .25 IN ZP 

i X/R -- "5.5IN 

1 UNIT HYPERBOLA -- .ujCH . . iWN 
ETAr. :3.8 







: — 1 


wm 



uTa 


z 



L . 

rxY^ 





rPK 

■mLw 


z 


rrxj\ 









rr 







O.L 0.2 0.-5 O.'u 0.5 o.b Q.7 O.M O.H i.O i,l .1.2 .1.3 i.i, .1,5 i.b 1.7 l.B i.^ 2 
f^REPAPfH BY S. KANYA AND D UIERENGA X/R 






APPENDIX 


1. OTHER POSSIBLE FORMS 


Wentzel (see References) cites the Schroedinger -Gordon wave 
equation, which has the form: 

(D (x,t) = 0, (1*1) 

0 = y 8^/9x^^ = (1.2) 

l/=l 


where 

V- I, 2, 3, 4, 

X , X , X = space coordinates, 

X M ^ 

X = i c t , 

4 


(1.3) 

(1.4) 

(1.5) 


it is known that the Broglie complex wave functions 

exp { i (k X + [p ^ + k^ c t (1.6) 

are particular solutions of the wave equation, (L 1)» The scalar con- 
stant, jLih c \ represents the rest mass of the respective particles. 

It turns out, indeed, that the stationary states of quantized fields 

obeying the field equations, (1. 1), represent systems of particles with 

-1 

the rest mass h c (Wentzel, p. 22). 

As is known, the association between quantum mechanics and 
macroscopic hydrodynamics in the classical sense was proposed by 
Irving Madelung in 1926. Exactly in the same year, Schroedinger 
announced his famous wave mechanics equation under the name of 
"amplitude equation." Madelung applies the "amplitude equation" of 
Schroedinger to the derivation of the equations of motion of the hydro - 
dynamic medium in the classical sense; the conditions and restrictions 
superimposed upon the Schroedinger equation and next superimposed 
by Madelung upon the hydrodynamic system in question and the manners 
of overcoming these difficulties are of the following nature: 

(a) The macroscopic hydrodynamic system in question when 
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moving in one direction should be free from the action of 
the curl. Of course one can use another system for the 
curl separately and add the results. 

(b) The Schroedinger equation refers to the amplitude and 
wave phenomena of one, single electron only. To use it 
in the sense of the macroscopic fluid dynamics, one has to 
propose some sort of generalization of the Schroedinger 
equation to a group or cluster of elements or molecules. 
Three such possible cases referring to the three possible 
groups of elements, h, m, were discussed in the present 
research. 

In that respect, the writer repeats the principles of the method from 

Methods of Quantum Field Theory in Statistical Physics , by A. A, 

Abrikosov, L. P. Gorkov, and I. E. Dzyalo shins ki, Institute for 

Physical Problems , Academy of Sciences, U.S.S.R., revised English 

edition, translated and edited by Richard A. Silverman; 

"In recent years, remarkable success has been achieved in 
statistical physics, due to the extensive use of methods 
borrowed from quantum field theory. The fruitfulness of 
these methods is associated with a new formulation of 
perturbation theory, primarily with the application of 
*Feynman diagrams. ' The basic advantages of the diagram 
technique lies in its intuitive character; Operating with 
one -particle concepts , we can use the technique to determine 
the structure of any approximation, and we can then write 
down the required expressions with the aid of correspondence 
rules. These new methods make it possible not only to solve 
a large number of problems which did not yield to the old 
formulation of the theory, but also to obtain many new re- 
lations of a general character. At present, these are the 
most powerful and effective methods available in quantum 
statistics . 

There now exists an extensive and very scattered journal 
literature devoted to the formulation of field theory methods 
in quantum statistics and their application to specific problems. 
However, familiarity with these methods is not widespread 
among scientists working in statistical physics. Therefore, 
in our opinion, the time has come to present a connected 
account of this subject, which is both sufficiently complete 
and accessible to the general reader." 

The present writer uses Feynman's {Nobel Prize) exact technique with 
the correspondence rules properly adjusted to the macroscopic fluid 
dynamics in two dimensions. The three-dimensional problems are not 
attacked, as yet, and have to be attacked in the future. 
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2. HEISENBERG OPERATORS 
2. 1. Preliminary Remarks 

There exists another formalism in attacking the problems in the 
field of quantum mechanics, i. e. , Heisenberg operators* This technique 
is not used for the time being in the present report and for this reason it 
will be only briefly mentioned here. May be in the future the writer may 
propose a method of application of Heisenberg operators to the problem 
which is under discussion in the present report. 

2, 2. Remarks on Operators 

Suppose that the function f{x) has a power series expansion: 

f(x) = f + X f , + x^f. + + . . . ; (2. 2. 1) 

' o I 2 3 

then, we can define the operator f{A) by: 

f(A) = f^ + A f ^ + A^f^ + A^f^ + . . . ; (2. 2. 2) 

for example, the operator exp (\A) is: 

exp (\A) = 1 + \ A + \^(2i )■ ^ A^ + )“ ^ A^ + . . . ; (2.2.3) 

it may be convenient now to introduce the proper notation used in the 
vector algebra: 

a column vector is denoted by: |*W> ; 

a row vector is denoted by; j ; 

examples: 

column vector: 

i’j'> = ^ I ; (2.2.4) 

a row vector: 

<4'| = .);*) , (2.2.5) 

where ^ stands for complex conjugate; the scalar product of a row 
vector [ and a column vector |^ > is equal to: 
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= <h iJj + d> lb — J (2, 2<p 6] 

y ^y 

if j\T'> has unit length then: 

|i|; + Itb 1^ = 1 ; (2.2.7; 

this is called the normalization condition^ As an example consider the 
following case: all beams of light are superpositions of many 

beams consisting of one photon each, one may turn his attention to the 
polarization properties of single photons. It is relatively easy to 
discover the probability rules for one photon from the knowledge of the 
behavior of classical beams. The general laws of quantum mechanics 
are just generalizations of these rules. For one photon one has {Baym, 
p. 3): 

lEp V = 8 TT te-u ; (2. 2. 8; 

where: 


E = E{r, t) = electric field vector; 

V = volume; 

4r = h(2iT) ^ ; 

w = the angular frequency. 

One defines the "state vector" of the photon polarization by: 


and by writing: 

= [V(8 ; 

■A. X 

= [V(87T-4rw) Ey ; 


the |^''> vectors are vectors in a "complex" two-dimensional space, 
their components being complex numbers. From Equation (2, 2, 7) it 
follows that > has unit length and: 

l+X = 1 . (2-2. 
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In fact the state vectors are independent of the volume V and depend only 
on the state of polarization of the photon. For example, if: 


^ > 


j / exp (i a ) 

^/2 \ exp (i a ) 


(2. 2. 13) 


then the photon is polarized at 45 to the x~axis, A knowledge of the 
[■^ > vector gives us all the information we can obtain about the state 
of the polarization of the photon. 

Some special examples of similar vectors are: 


1 X> 

^ 0 

: x-polarization; 

(2.2.14) 

y > 

= 0 

: y- polarization ; 

(2. 2. 15) 

R > 

1_ 

^/2 

[ ^ ; right circular polarization; 

/ 

(2. 2. 16) 

1 L> 

1 

nT2 

^ j ; left circular polarization . 

(2.2. 17) 


Let us associate with each column vector | > a row vector < "4^ \ 

which is defined by: 


where stands for complex conjugate. We may define some operations 
like the scalar product (discussed already above): scalar product of a 

row vector. 1 and a column vector | ;^ > : 


< ^ I ^ > 


<|> + 


<j) Jll 

y 


= < 


^ 1 $>‘ 


(2.2. 19) 


the known normalization condition requires for the vector 




I ^ > = 1 ; 

clearly: 

<x |x> =<y|y> = 1; 


( 2 . 2 . 20 ) 


( 2 . 2 . 21 ) 
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<r| 

R> = <L 

V 

II 

(2, 2. 22) 

the vectors 

|x > 

and 1 y > 

are orthogonal, i. e. 

, perpendicular, in the 

sense that: 

<x| 

o 

It 

A 


(2, 2. 23) 

The vectors 

1 x> 

and 1 y > 

form a basis, since 

any vector can be 


written as a linear superposition of them: 


^ > 


/ v' 

\*yl 


= I ^ ^ 1 y > 


(2. 2o 24) 


they are orthogonal, they satisfy the normalization condition, and they 
form the ' 'orthonormal" basis. Similarly, the set | R> and |L* > 
form an orthonormal basis, since we can always write: 


^ > 


IK\ 

\^l 


X X |r> + XE Y- |l> 


sfZ 




(2. 2. 25) 


We can regard any arbitrary polarization as a coherent (logically 
consistent) superposition of x and y polarization states, or equivalently 
as a coherent superposition of right and left circularly polarization states. 
Let us return to the problem of passing a beam through an x-polaroid. 

The classical rules tell us to regard the beam as a superposition of an 
x-polarized beam and a y-polarized beam, and that the effect of a 
Polaroid is to throw away the y-polarized component and pass only the 
x-polarized component. The absolute value squared of the amplitude 
of the beam gives us its energy before it passes through, and the 
absolute value squared of its x- component give us its energy after it 
passes through. The fraction of the beam that passes through is given 
by equation: 


1 = |E^|^ [ |E|^r' ; (2.2.26) 


as we know, quantum mechanically the fraction (2, 2. 26) gives us the 
probability of one photon with the initial polarization passing through the 
Polaroid, Written in terms of |»>|j > , Eq, (2, 2. 26) is the "probability": 
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probability = |i]j |^ [ ^ ” I I • 2*27) 

Thus < X > is the amplitude of the x- polarized component of |'T/> , 
and its absolute value squared is the probability that the photon in the 
state > passes through the x~polaroid* We call <xj"5^ > the 
’’probability amplitude" for the photon to pass through the x-polarizer. 

Next, consider passing light through a prism that passes right 
circular light but rejects left circular light. To calculate the probability, 
we write the beam as a coherent sum of right circularly polarized and 
left circularly polarized light, in terms of the state vectors. Then 
<r|'^> is the amplitude of the component passed through the prism and 
I <R 1^ is the probability that a photon in the state > will pass 
through the prism that passes only the right circular polarized light. 

Thus the general rule is that if we have a prism that passes only light 
in the state > , rejecting light in states orthogonal to |^> , then the 

probability amplitude that a photon in the state will pass through 

the prism is; 

<$ 1^ > , (2. 2. 28) 
and the probability that the photon passes through is: 

|<$|^>|2 . (2.2.29) 

Notice that the probability is independent of the phase of |^ > or [0> . 

2. 3. Equations 

Above, we defined the operator f(A) as: 

f{A) = f^ + A f j + A^f^ + A^f^ + ... . (2.3.1) 

An application of an operator (2. 3. 1) to any "state vector," Equation 
(2. 2. 9), say, causes a transformation of this vector. The operator, 
composed of the right and left circular polarizations, is usually denoted 
by the symbol R(0) . In general, if we transform a vector |^> by R(0), 
the new components of |^ > will bear little resemblance to the old 
components. Let us ask if there are any state vectors that when 
operated upon by R(0) are at most multiplied by a constant; 
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R(0)|^ > = C |vp > . 


(2. 3. 2) 


If a vector satisfies a relation like this it is called an ''eigenvector*' of 
R(0), and the number C is called the "eigenvalue” of R{9) belonging to 
the eigenvector in questiono Analogously, one can talk about the "eigen- 
states” of an operator* If |a> is an eigenstate of f(A), Equation (2, 3* 1), 
with eigenvalue a, then operating on |a> with f(A), one gets: 

f (A) la > = (f^ + A f ^ d A^f^ + A^f^ + * o . ) |a > 

= (f^ + a f j o » • ) l<i> 

= f(a) la > o (2* 3o 3 ) 


The above operation emphasizes clearly that |a> is an eigenstate of the 
operator f(A) with an eigenvalue f(a)o The eigenstates of an operator A 
form the so-called "complete” set of eigenstates of A; one of the conditions, 
to be satisfied by the set in question is the following one: 

l = ^la><a|; (2. 3. 4) 

a 


the eigenstates of A can always be chosen so that a complete set of A, 
they form is Hermitian* Operating on both sides of Equation (2* 3* 4) with 
f(A) one finds from Equation (2. 3* 4); 


f(A) 


^ f(a)(a><a 
a 


(2, 3* 5) 


In fact, one could take this (i. e* , Equation (2* 3» 5)) as the equation 
defining f{A) in terms of the eigenstates of A; it is not necessary for ”f” 
to have a power series expansion to be a usable definition* For illustrative 
purposes, one can define a delta function of an operator, 6(A-x), by: 


6(A-x) 




6(a-x) a> < a 


( 2 . 3 . 6 ) 


where x is number (Baym, p. 13 5). 

To complete the explanations of the notation, the writer puts down 
a few more notations used in the text below. Let the letter "i” denote 
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the length in the z direction, say, on which length some changes in the 
states of a phenomenon may occur. Let the incident state in a phenomenon 
in question be denoted by the symbol given below: 

incident state = je > ; (2. 3. 7) 

the distance i changes the incident state |e> into the state: 

je > exp (i k^i ) |e > ; (2. 3. 8) 

where: 

k = wave number; ('^ cm ^); 
e 

example in the electric field: 

E (r*, t) = cos (k z “ wt + a ) ; (2. 3. 9) 

XX X 


E (r^,t) = cos (k z - ut + a ) ; (2. 3. 10) 

y' ' ^ y ' y' 

k = 2tt X~ ^ = wave number; '^(cm ^); 
u = angular frequency ; (^ sec ^); 

period of a vibratory motion = time required for a complete vibration; 
frequency = number of complete vibrations, or cycles, in a unit of time; 
the frequency is numerically equal to the reciprocal of the period; 

X. = wave length of the harmonic wave = distance between successive 
points at which the wave differs in phase by 2tt, or the distance between 
matximum displacements or crests; a , a = the phases; phase is a 
term stating where the body is at any instant in its vibrations, and what 
its direction of motion is; the wave numbers with subscripts denote: 


, -1 / - 1 ^ 

k = w n c ; (~ cm ); 

e e 

k = w n c ; (~ cm ); 

o o 


(2. 3. 11) 
(2. 3. 12) 


where n , n , are integers. The wave numbers often do not have any 
e o 

special names. 

If the phenomenon in question takes place on the length *'i", then 
the difference in phase of an ”e wave" between its point of entry and its 
point of departure will be "k i Similarly, an "O wave" changes phase 
by "k^i". The time factor exp(- iiot) is the same for both states (of 
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polarization, say) and does not effect their relative phase. In Equation 
(2. 3. 8) the incident state le> is changed on the distance jl ; this implies 
tliat (as mentioned already in Equation (2, 3. 8))t 


e > exp {i k^i ) je > ; 


(2, 3, 13) 


as it is seen, the function, which introduces (causes) the changes in the 
states is assumed to be in the form, of the exponential function, exp (i k^i ), 
Similarly, and incident jo^ state changes into 

jo > exp (i k^i ) |0 > ; (2. 3, 14) 

the state jO > is associated with the wave number k^, i, e. , : 

state Jo > ; wave number k^ = w ^ ; (2, 3, 15) 

let us define the wave number "matrix" by: 


K = k |e> <e| + k„|0x0l , (2,3.16) 

e u 

-f- 

The Hermitian adjoint, K , of the matrix K, Equation (2, 3, 16), is the 
matrix: 



1 <x|k|x>’ 

II * ^ 

<y )K|x> 

K+ = 


»JO 


^<x|K|y>‘' 

<y|K|y>" 

/ 

It may happen that in some cases (see Baym 


(2, 3, 17) 


K' 


= K 


(2. 3, 18) 


K is then said to be Hermitian, 

As an application of operator functions, one may cite the following 
example and result: 

Suppose that |^(0)> is the state of a system at time t = 0; then 
at another time t, the state of the system is given by the equation: 



exp {- i Et^“^lE><El^(0)> 


(2. 3. 19) 


where the sum is over a complete set of energy eigenstates of the Hamiltonian, 
H, of the system. 
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The row vector <S^(t)| is given in terms of <^(0) | by: 

< ^(t) I = <'lf{0)| (exp(iHt-ii 

= <^(0)| exp (i Ht -hr , (Z. 3. 20) 

where the hermiticity of H is used (see Baym, p, 136). Using Equation 
(2. 3. 5), Equation (2. 3. 19) is simply: 

|^(t) > = exp (“ i Ht4r |'®^(0)> ; (2. 3. 21) 

equations (2. 3. 19) and (2. 3. 21) are only valid if H does not explicitly 
depend on time (Baym, p, 135). The expression |‘^(t)> in this form 
clearly solves the Schroedinger equation: 

i h =Hl^(t)> . (2.3.22) 

dt 

2. 4. Heisenberg Representation 

On may be often interested in knowing how expectation values of 
operators, in the state |^(t)> , change in course of time. If A is an 
operator, then its expectation value at time t is: 

<A>^ = <^(t)|A|^(t) > . (2.4.1) 

Assume A not to depend on time explicitly. Using Equation (2. 3. 20) and 
(2. 3. 21) one gets: 

< A>^ = <^(0) |exp (i Htlr ^) A exp (+i Htlr ^ | ^(0) > ; (2.4.2) 

if we define a time- dependent operator A(t) by: 

A(t) = exp (i Ht 4r ^) A exp (- i Ht4r ^) , (2, 4, 3) 

then we can also write; 

<A>^ = <^(0)|A(t)|^(0) > . (2.4.4) 

In this equation one can regard the time development of <A>^ as occurring 
because the operator changes in time, while the state remains the same at all 
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times. This way of regarding the time development of the system is called 
the ''Heisenberg representation" or "Heisenberg picture." The other 
v/ay of looking at the time development, as in Equation .C Z, 3 , 22 ), where 
the operators remain constant in time, but the states change according to 
Equation ( 2.3.22 ), is called the "Schroedinger representation" (or 
picture). At t = 0, the states and operators are the same in both repre- 
sentations. Both representations give the same results for time dependent 
expectation values. We can solve for either the time dependence of the 
states in the Schroedinger representation, or for the time dependence of 
the operators in the Heisenberg representation. 

In the present report the writer applies the Schroedinger represen- 
tation to the problem proposed and outlined by Madelung, i.e. , the 
association between the quantum mechanics and macroscopic hydro- 
dynamics . 
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